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SUMMARY 


This  report  documents  the  application  of  a  multigrid  Computational  Fluid  Dynamics 
(CFD)  method  to  rotor  aerodynamic  prediction.  The  CFD  code  TLNS3DR  is  used.  The 
code  is  the  modified  version  of  the  nonrotating  wing  TLNS3D  code  with  the  addition  of 
the  rotation  term.  The  code  solves  the  three-dimensional  thin-layer  Navier-Stokes  equations 
using  an  explicit  multistage  Runge-Kutta  type  of  time  stepping  with  multigrid  method.  For 
nonrotating  wings,  the  TLNS3DR  code  is  identical  to  the  TLNS3D  code. 

The  TLNS3DR  code  is  applied  to  various  advanced  rotor  blades  to  study  tip  vortex 
and  to  investigate  the  capability  of  the  code  in  rotor  analysis.  For  nonrotating  flows,  the 
solutions  in  terms  of  the  tip  vortex  and  surface  pressure  coefficients  are  obtained  for  realistic 
helicopter  rotor-tip  configurations  to  study  effects  of  blade  planforms  on  the  tip  vortex  under 
incompressible  flow  condition.  Calculated  results  are  compared  with  experimentally  obtained 
data  at  NASA-Langley’s  Basic  Aerodynamic  Research  Wind  Tunnel  (BART).  The  Berp-tip, 
the  swept-tip  and  an  equivalent  taper-tip  configurations  are  investigated.  For  both  nonrotating 
and  rotating  flows,  the  Berp-type  blade,  the  swept-type  blade,  and  equivalent  taper-type 
and  rectangle-type  blades  with  zero  twist  are  considered  under  compressible  flow  condition. 
Solutions  for  rotating  flows  are  presented  along  with  those  of  nonrotating  flows  to  study  the 
effect  of  the  rotation  term  added  to  the  TLNS3DR  code  in  addition  to  the  study  of  the  tip 
vortex. 


CONTENTS 


Page 

PREFACE  .  iii 

SUMMARY  .  v 

LIST  OF  FIGURES  .  viii 

1.  INTRODUCTION  .  1 

2.  METHODOLOGY  OF  THE  TLNS3DR  CODE  . 2 

2.1  NS  Equations  in  Blade-fixed  Frame  .  2 

2.2  Finite- Volume  Discretization  .  4 

2.3  Multi-stage  Time-Stepping  Scheme  .  4 

2.4  Multigrid  Acceleration  Technique  .  5 

2.5  Boundary  Conditions  and  Grid  .  5 

3.  COMPUTATIONAL  RESULTS  .  6 

3.1  Nonrotating  Blade  Tips  .  6 

3.1.1  Tip-flow  field  at  a  =  15°  6 

3.1.2  Tip-flow  field  at  a  =  5°  —  7.5°  7 

3.1.3  Tip-flow  field  at  a  =  0°  7 

3.1.4  Surface  pressure  and  convergence  history  .  7 

3.2  Comparison  of  Rotating  and  Nonrotating  Rotor  Flows  .  8 

4.  CONCLUDING  REMARKS  .  9 

REFERENCES  .  144 

SYMBOLS  .  145 


FIGURES 


Figure  Page 

1  Plane  form  of  Berp-tip  .  10 

2  Plane  form  of  swept-tip  . .  H 

3  Plane  form  of  taper-tip  .  . .  12 

4  View  of  the  twist  for  Berp-tip  .  13 

5  View  of  the  twist  for  swept-tip  .  14 

6  View  of  the  twist  for  taper-tip  .  15 

7  Partial  view  of  C-0  mesh  for  Berp-tip  .  16 

8  Partial  view  of  C-0  mesh  for  swept-tip  .  17 

9  Partial  view  of  C-0  mesh  for  taper-tip  .  18 

10  Calculated  tip  flow  field  in  terms  of  downwash  for  Berp-tip  at  a  =  15°  19 

11  Calculated  tip  flow  field  in  terms  of  downwash  for  swept-tip  at  a  =  15°  20 

12  Calculated  tip  flow  field  in  terms  of  downwash  for  taper-tip  at  a  =  15°  21 

13  Calculated  tip  flow  field  in  terms  of  downwash  for  taper-tip  at  a  =  15° 

with  a  fine  grid  .  22 

14  Measured  tip  flow  field  in  terms  of  downwash  for  Berp-tip  at  a  =  15°  .  23 

15  Measured  tip  flow  field  in  terms  of  downwash  for  swept-tip  at  a  =  15°  24 

16  Measured  tip  flow  field  in  terms  of  downwash  for  taper-tip  at  a  =  15°  25 

17  Calculated  tip  flow  field  in  terms  of  axial  velocity  for  Berp-tip  at  a  =  15°  •  •  •  •  26 

18  Calculated  tip  flow  field  in  terms  of  axial  velocity  for  swept-tip  at  a  =  15°  •  •  •  27 

19  Calculated  tip  flow  field  in  terms  of  axial  velocity  for  taper-tip  at  a  =  15°  •  •  •  28 

20  Calculated  tip  flow  field  in  terms  of  axial  velocity  for  taper-tip  at  a  =  15° 

with  a  fine  grid  .  29 

21  Measured  tip  flow  field  in  terms  of  axial  velocity  for  Berp-tip  at  a  =  15°  .  30 

22  Measured  tip  flow  field  in  terms  of  axial  velocity  for  swept-tip  at  a  =  15°  •  •  •  •  31 

23  Measured  tip  flow  field  in  terms  of  axial  velocity  for  taper-tip  at  a  =  15°  ....  32 

24  Calculated  tip  flow  field  in  terms  of  crossflow  velocity  vector  for  Berp-tip 


vm 


at  a  =  15°  . 

25  Calculated  tip  flow  field  in  terms  of  crossflow  velocity  vector  for  swept-tip 

at  a  =  15°  . 

26  Calculated  tip  flow  field  in  terms  of  crossflow  velocity  vector  for  taper-tip 

at  a  =  15°  . 

27  Calculated  tip  flow  field  in  terms  of  crossflow  velocity  vector  for  taper-tip 

at  a  =  15°  with  a  fine  grid  . 

28  Measured  tip  flow  field  in  terms  of  crossflow  velocity  vector  for  Berp-tip 

at  a  =  15°  . 

29  Measured  tip  flow  field  in  terms  of  crossflow  velocity  vector  for  swept-tip 

at  a  =  15°  . 

30  Measured  tip  flow  field  in  terms  of  crossflow  velocity  vector  for  taper-tip 

at  a  =  15°  . 

31  Calculated  tip  flow  field  in  terms  of  downwash  for  Berp-tip  at  a  =  5°  . 

32  Calculated  tip  flow  field  in  terms  of  downwash  for  swept-tip  at  a  =  5°  . 

33  Calculated  tip  flow  field  in  terms  of  downwash  for  taper-tip  at  a  =  7.5°  . 

34  Calculated  tip  flow  field  in  terms  of  axial  velocity  for  Berp-tip  at  a  =  5°  .  .  .  . 

35  Calculated  tip  flow  field  in  terms  of  axial  velocity  for  swept-tip  at  a  =  5°  •  •  • 

36  Calculated  tip  flow  field  in  terms  of  axial  velocity  for  taper-tip  at  a  =  7.5°  •  •  • 

37  Calculated  tip  flow  field  in  terms  of  crossflow  velocity  vector  for  Berp-tip 

at  a  =  5°  . 

38  Calculated  tip  flow  field  in  terms  of  crossflow  velocity  vector  for  swept-tip 

at  a  =  5°  . 

39  Calculated  tip  flow  field  in  terms  of  crossflow  velocity  vector  for  taper-tip 

at  a  =  7.5°  . 

40  Measured  tip  flow  field  in  terms  of  crossflow  velocity  vector  for  three  tips 

at  a  =  5°  to  7.5°  and  x/c  =  1.65  . 

41  Calculated  tip  flow  field  in  terms  of  downwash  for  Berp-tip  at  a  =  0°  . 

42  Calculated  tip  flow  field  in  terms  of  downwash  for  swept-tip  at  a  =  0°  . 

43  Calculated  tip  flow  field  in  terms  of  downwash  for  taper-tip  at  a  =  0°  . 

44  Calculated  tip  flow  field  in  terms  of  downwash  for  taper-tip  at  a  =  0° 

with  a  fine  grid  . 

45  Calculated  tip  flow  field  in  terms  of  axial  velocity  for  Berp-tip  at  a  =  0°  •  •  •  • 


33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 


IX 


55 


46  Calculated  tip  flow  field  in  terms  of  axial  velocity  for  swept-tip  at  a?  =  0° 

47  Calculated  tip  flow  field  in  terms  of  axial  velocity  for  taper-tip  at  a  =  0°  ....  56 

48  Calculated  tip  flow  field  in  terms  of  axial  velocity  for  taper-tip  at  a  =  0° 

with  a  fine  grid  .  57 

49  Calculated  tip  flow  field  in  terms  of  crossflow  velocity  vector  for  Berp-tip 

at  a  =  0°  .  58 

50  Calculated  tip  flow  field  in  terms  of  crossflow  velocity  vector  for  swept-tip 

at  of  =  0°  .  59 

51  Calculated  tip  flow  field  in  terms  of  crossflow  velocity  vector  for  taper- tip 

at  a  =  0°  .  60 

52  Calculated  tip  flow  field  in  terms  of  crossflow  velocity  vector  for  taper-tip 

at  a  =  0°  with  a  fine  grid  .  61 

53  Measured  tip  flow  field  in  terms  of  crossflow  velocity  vector  for  three  tips 

at  a  =  0°  and  x/c  —  1.65  .  62 

54  Calculated  surface  pressure  coefficients  for  Berp-tip  at  a  =  15°  63 

55  Calculated  surface  pressure  coefficients  for  swept-tip  at  a  =  15°  64 

56  Calculated  surface  pressure  coefficients  for  taper-tip  at  a  =  15°  65 

57  Calculated  surface  pressure  coefficients  for  Berp-tip  at  a  =  5°  66 

58  Calculated  surface  pressure  coefficients  for  swept-tip  at  a  =  5°  .  67 

59  Calculated  surface  pressure  coefficients  for  taper-tip  at  a  =  7.5°  68 

60  Calculated  surface  pressure  coefficients  for  Berp-tip  at  a  =  0°  69 

61  Calculated  surface  pressure  coefficients  for  swept-tip  at  a  =  0°  70 

62  Calculated  surface  pressure  coefficients  for  taper-tip  at  a  =  0°  71 

63  Convergence  history  in  terms  of  residuals  for  Berp-tip 

at  a  =  15°,  a  =  5°  and  a  =  0°  .  72 

64  Convergence  history  in  terms  of  residuals  for  swept-tip 

at  a  =  15°,  a  —  5°  and  a  =  0°  .  73 

65  Convergence  history  in  terms  of  residuals  for  taper-tip 

at  a  —  15°,  a  =  7.5°  and  a  =  0°  .  74 

66  Convergence  history  in  terms  of  residuals  for  taper-tip 

at  a  =  15°  and  a  =  0°  with  a  fine  grid  . .  75 

67  Plane  form  of  Berp-blade  .  76 


x 


68  Partial  view  of  C-0  mesh  for  Berp-blade  .  77 

69  Plane  form  of  swept -blade  .  78 

70  Partial  view  of  C-0  mesh  for  swept-blade  .  79 

71  Plane  form  of  taper-blade  .  80 

72  Partial  view  of  C-0  mesh  for  taper-blade  .  81 

73  Plane  form  of  rectangle-blade  . . .  82 

74  Partial  view  of  C-0  mesh  for  rectangle-blade  .  83 

75  Calculated  tip  flow  field  in  terms  of  axial  velocity  for 

nonrotating  Berp-blade  .  84 

76  Calculated  tip  flow  field  in  terms  of  crossflow  velocity  vector  for 

nonrotating  Berp-blade  .  85 

77  Calculated  surface  pressure  coefficients  for  nonrotating  Berp-blade  .  86 

78  Calculated  local  Mach  contours  for  nonrotating  Berp-blade  .  87 

79  Calculated  local  V-velocity  contours  for  nonrotating  Berp-blade  .  88 

80  Convergence  history  in  terms  of  residuals  for  nonrotating  Berp-blade  .  89 

81  Calculated  tip  flow  field  in  terms  of  relative  axial  velocity  for 

Berp-blade  in  hover  .  90 

82  Calculated  tip  flow  field  in  terms  of  absolute  axial  velocity  for 

Berp-blade  in  hover  .  91 

83  Calculated  tip  flow  field  in  terms  of  relative  crossflow  velocity  vector  for 

Berp-blade  in  hover  .  92 

84  Calculated  tip  flow  field  in  terms  of  absolute  crossflow  velocity  vector  for 

Berp-blade  in  hover  .  93 

85  Calculated  surface  pressure  coefficients  for  Berp-blade  in  hover  .  94 

86  Calculated  local  relative  Mach  contours  for  Berp-blade  in  hover  .  95 

87  Calculated  local  relative  V-velocity  contours  for  Berp-blade  in  hover  .  96 

88  Calculated  local  absolute  V-velocity  contours  for  Berp-blade  in  hover  .  97 

89  Convergence  history  in  terms  of  residuals  for  Berp-blade  in  hover  .  98 

90  Calculated  tip  flow  field  in  terms  of  axial  velocity  for 

nonrotating  swept-blade  . . .  99 


xi 


91  Calculated  tip  flow  field  in  terms  of  crossflow  velocity  vector  for 

nonrotating  swept-blade  .  100 

92  Calculated  surface  pressure  coefficients  for  nonrotating  swept-blade  .  101 

93  Calculated  local  Mach  contours  for  nonrotating  swept-blade  .  102 

94  Calculated  local  V-velocity  contours  for  nonrotating  swept-blade  .  103 

95  Convergence  history  in  terms  of  residuals  for  nonrotating  swept-blade  .  104 

96  Calculated  tip  flow  field  in  terms  of  relative  axial  velocity  for 

swept-blade  in  hover  .  105 

97  Calculated  tip  flow  field  in  terms  of  absolute  axial  velocity  for 

swept-blade  in  hover  .  106 

98  Calculated  tip  flow  field  in  terms  of  relative  crossflow  velocity  vector  for 

swept-blade  in  hover  .  107 

99  Calculated  tip  flow  field  in  terms  of  absolute  crossflow  velocity  vector  for 

swept-blade  in  hover  .  108 

100  Calculated  surface  pressure  coefficients  for  swept-blade  in  hover  .  109 

101  Calculated  local  relative  Mach  contours  for  swept-blade  in  hover  .  110 

102  Calculated  local  relative  V-velocity  contours  for  swept-blade  in  hover  .  Ill 

103  Calculated  local  absolute  V-velocity  contours  for  swept-blade  in  hover  .  112 

104  Convergence  history  in  terms  of  residuals  for  swept-blade  in  hover  .  113 

105  Calculated  tip  flow  field  in  terms  of  axial  velocity  for 

nonrotating  taper-blade  . . .  114 

106  Calculated  tip  flow  field  in  terms  of  crossflow  velocity  vector  for 

nonrotating  taper-blade  .  115 

107  Calculated  surface  pressure  coefficients  for  nonrotating  taper-blade  .  116 

108  Calculated  local  Mach  contours  for  nonrotating  taper-blade  .  117 

109  Calculated  local  V-velocity  contours  for  nonrotating  taper-blade  .  118 

110  Convergence  history  in  terms  of  residuals  for  nonrotating  taper-blade  .  119 

111  Calculated  tip  flow  field  in  terms  of  relative  axial  velocity  for 

taper-blade  in  hover  .  120 

112  Calculated  tip  flow  field  in  terms  of  absolute  axial  velocity  for 

taper-blade  in  hover  .  121 

113  Calculated  tip  flow  field  in  terms  of  relative  crossflow  velocity  vector  for 

xii 


taper-blade  in  hover  .  122 

114  Calculated  tip  flow  field  in  terms  of  absolute  crossflow  velocity  vector  for 

taper-blade  in  hover  . 123 

115  Calculated  surface  pressure  coefficients  for  taper-blade  in  hover  .  124 

116  Calculated  local  relative  Mach  contours  for  taper-blade  in  hover  .  125 

117  Calculated  local  relative  V-velocity  contours  for  taper-blade  in  hover  .  126 

118  Calculated  local  absolute  V-velocity  contours  for  taper-blade  in  hover  .  127 

119  Convergence  history  in  terms  of  residuals  for  taper- blade  in  hover  .  128 

120  Calculated  tip  flow  field  in  terms  of  axial  velocity  for 

nonrotating  rectangle-blade  .  129 

121  Calculated  tip  flow  field  in  terms  of  crossflow  velocity  vector  for 

nonrotating  rectangle-blade  . 130 

122  Calculated  surface  pressure  coefficients  for  nonrotating  rectangle-blade  ....  131 

123  Calculated  local  Mach  contours  for  nonrotating  rectangle-blade  .  132 

124  Calculated  local  V-velocity  contours  for  nonrotating  rectangle-blade  .  133 

125  Convergence  history  in  terms  of  residuals  for  nonrotating  rectangle-blade  •  •  •  134 

126  Calculated  tip  flow  field  in  terms  of  relative  axial  velocity  for 

rectangle-blade  in  hover  •  •  . .  135 

127  Calculated  tip  flow  field  in  terms  of  absolute  axial  velocity  for 

rectangle-blade  in  hover  .  136 

128  Calculated  tip  flow  field  in  terms  of  relative  crossflow  velocity  vector  for 

rectangle-blade  in  hover  .  137 

129  Calculated  tip  flow  field  in  terms  of  absolute  crossflow  velocity  vector  for 

rectangle- blade  in  hover  .  138 

130  Calculated  surface  pressure  coefficients  for  rectangle-blade  in  hover  .  139 

131  Calculated  local  relative  Mach  contours  for  rectangle-blade  in  hover  .  140 

132  Calculated  local  relative  V-velocity  contours  for  rectangle-blade  in  hover  .  .  .  141 

133  Calculated  local  absolute  V-velocity  contours  for  rectangle-blade  in  hover  .  .  142 

134  Convergence  history  in  terms  of  residuals  for  rectangle-blade  in  hover  .  143 

xiii 


1.  INTRODUCTION 


The  accurate  and  efficient  prediction  of  rotor  flow  fields  poses  a  particularly  challenging 
problem  in  helicopter  rotor  design.  The  accurate  prediction  of  the  aerodynamic  forces  on  the 
rotor  requires  a  correct  prediction  of  the  vortex  wake,  including  the  blade  trailing  wake  as  well 
as  the  blade- vortex  interaction  (BVI).  The  BVI  noise  is  also  one  of  the  most  important  noise 
sources  for  rotary  wings.  The  detailed  experimental  test  for  predicting  rotor  flow  field  is  very 
expensive  and  complex,  which  limits  the  amount  of  available  data,  particularly  on  advanced 
rotors.  Therefore,  other  methods  such  as  Computational  Fluid  Dynamics  (CFD)  are  needed 
in  rotor  analysis. 

The  CFD  methods  for  calculating  helicopter  rotor  flow  fields  can  be  divided  into  two 
types:  integral  equation  methods  (panel  methods)  and  finite-difference/finite- volume  methods. 
Integral  equation  methods  are  restricted  to  linear,  low-speed,  subsonic,  or  at  most,  to  nonlinear 
flows  with  weak  shocks.  This  is  a  very  serious  limitation,  because  modern  rotor  blades  operate 
in  the  transonic  regime,  where  the  nonlinear  compressibility  is  important  and  shocks  can 
exist  in  the  flow  field.  Finite-difference/finite- volume  methods  based  on  Navier-Stokes  (NS) 
equations  with  a  single  grid  for  rotor  flow  fields  were  developed1,2.  The  methods  can  treat 
compressible  flows  with  shocks,  including  strong  shocks,  in  the  region  of  the  rotor  blade 
accurately.  However,  the  calculation  of  NS  equations  on  a  fine,  single  grid  is  expensive. 

A  multigrid  NS  CFD  computer  code,  TLNS3D,  was  developed3  for  the  fixed-wing  motion. 
The  TLNS3D  code  solves  the  three-dimensional  thin-layer  NS  equations  using  an  explicit 
multistage  Runge-Kutta  type  of  time-stepping  with  multigrid  method.  The  philosophy  of 
the  multigrid  method  is  the  use  of  successively  coarser  grids  to  compute  corrections  to  a 
fine  grid  solution.  The  use  of  multigrid  strategy  results  in  a  much  faster  convergence  than 
the  single-grid  method.  The  efficiency  of  the  multigrid  method  using  five-stage  Runge-Kutta 
time-stepping  was  demonstrated3  through  computations  for  transonic  flows  over  fixed  wings. 
The  code  was  shown  to  result  in  a  significant  reduction  in  CPU  time  as  compared  to  the 
single-grid  multistage  time-stepping  method.  The  code  is  recently  modified  for  rotary-wing 
calculations4  with  the  addition  of  the  rotation  term,  although  the  accuracy  of  the  code  for 
rotary- wing  calculations  has  not  been  tested.  The  new  code  called  TLNS3DR  is  identical  to 
the  TLNS3D  code  for  nonrotating  wings. 

The  purpose  of  the  present  work  is  to  use  the  efficient  multigrid  TLNS3DR  CFD  code 
to  analyze  advanced  rotor  tip  designs  and  to  further  investigate  the  capability  of  the  code  for 
rotor  aerodynamic  analysis.  For  nonrotating  flows,  the  solutions  in  terms  of  the  tip  vortex 
and  surface  pressure  coefficients  are  obtained  for  realistic  helicopter  rotor-tip  configurations 
to  study  the  effects  of  blade  planforms  on  the  tip  vortex  under  incompressible  flow  condition. 
Calculated  results  are  compared  with  experimentally  obtained  data  at  NASA-Langley’s  Basic 
Aerodynamic  Research  Wind  Tunnel  (BART).  The  Berp-tip,  the  swept-tip  and  an  equivalent 
taper-tip  configurations  are  investigated.  For  both  nonrotating  and  rotating  flows,  the  Berp- 
type  blade,  the  swept-type  blade,  and  equivalent  taper-type  and  rectangle-type  blades  with 
zero  twist  are  considered  under  compressible  flow  condition.  Solutions  for  rotating  flows  are 
presented  along  with  those  of  nonrotating  flows  to  study  the  effect  of  the  rotation  term  added 
to  the  TLNS3DR  code  in  addition  to  the  tip  vortex  study. 
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2.  METHODOLOGY  OF  THE  TLNS3DR  CODE 


2.1  NS  Equations  in  Blade-fixed  Frame 

For  a  general  motion  of  a  rotor  blade,  the  governing  equations  axe  simple  to  solve  if 
a  rotor  blade-fixed  moving  frame  of  reference  formulation  is  used.  The  moving  frame  of 

reference  (x,  y,  z )  is  translating  at  a  velocity  of  V0(t)  and  rotating  around  a  pivot  point  (rotor 
blade  axis^  rp  =  (xp,yp,  zp))  at  an  angular  velocity  of  Q,(t).  The  relation  for  the  absolute 
velocity  (V  =  (u,  relative  velocity  (Vr  =  («r,  ty,  wrj)  and  transformation  velocity 

(Vt  =  (uf ,  wt )  =  V0  +  Ve  =  Va  +  Q,  x  r)  is  given  by 

V  =  Vr  +  V0  +  Q  x  r  (1) 

where  r  is  the  position  vector  measured  in  the  blade-fixed  moving  frame  of  reference  relative 
to  Tp.  It  should  be  noted  that  the  transformation  velocity  (Vt)  includes  both  translation  (V0) 

and  rotation  (fi  x  r)  velocities  and  therefore  this  formulation  is  general  for  both  hovering  and 
forwarding  motion  of  the  helicopter.  For  the  fixed-wing  motion,  the  transformation  velocity 
is  set  to  zero,  and  hence  the  relative  velocity  is  the  same  as  the  absolute  velocity. 

The  NS  equations  in  terms  of  the  relative  velocity  in  the  blade-fixed  moving  frame  of 
reference  is  specialized  to  a  body  fitted  coordinate  system  (£,77,  £)  in  computational  space, 
where  f,  77  and  (  represent  the  streamwise,  normal  and  spanwise  directions,  respectively. 
The  thin-layer  assumption  is  employed  by  keeping  a  viscous  diffusion  term  in  the  normal  (77) 
direction  only,  since  the  dominant  viscous  effects  arise  from  viscous  diffusion  normal  to  the 
body  surface  for  high-Reynolds-number  turbulent  flows.  For  the  time-independent  (£,  77,  () 
coordinate  system  in  the  blade-fixed  moving  frame  of  reference,  the  thin-layer  NS  equations 
take  the  following  form3,4: 


with 


&('-1Vr)  +  %  +  ^+dHr 


d£  drj  '  d( 


dGm 

drj 


+  J~1S 


Ur  =  [p,pUr,pVr,pWr,pe  r]* 

Fr  =  J~l  [pur,  purur  +  (xp,  purVr  +  tyP,  purwr  +  (zp,  purhry 
Gr  =  J~1  [pVr,  pVrUr  +  77 x p ,  pVrVr  +  T)yp,  pVrWr  +  T]zP,  P^rK]* 
Hr  —  J  1  [pWr,  pU)rUr  +  £ xP ,  pWrVr  +  ( yp ,  pWrWr  +  (zp,  pWrhrY 


_  j-xyfyMoope 

=  J  ft  [0)  </>lurT,  +  Th<t>2,<f>lVrT,  +  Tlyfa,  <t>\wrv  +  Vz<t>2i  <f>l  a  +  Vrfc]* 


S  =  [0,  -patx ,  -paty ,  -patz , 

-  p[Vr  ■a0  +  (Q,xr)-a0  +  V0-(at-&xVr) 

+  '  (jF  x  ^  +  ^  x  ^ x 

where 

<h  =  vl  +  vl  +  T)] 


(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 


(9) 
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(10) 


<h  =  +  riyV r,  +  l),Wr,) 

<n> 

?2  =  it2  +  v\  +  wl  (12) 

where  the  prime  (/)  denotes  the  operation  with  respect  to  the  blade-fixed  moving  frame  of 
reference.  In  Eqs.  (7)  and  (11),  the  turbulent  eddy-viscosity  correction  e  and  turbulent  eddy- 
conductivity  correction  e  are  used.  The  algebraic  turbulent  model  of  Baldwin  and  Lomax5  is 
employed. 


The  pressure  p  and  the  total  enthalpy  hr  are  calculated  by 

V2  V2 

p  =  (7  -  lWr  -  -i-  +  -f) 


(13) 


and 

h  _  7 P  +Vr2  V2 
"~p(  7-l)+2  2 


(14) 


respectively.  The  at  is  the  acceleration  of  the  transformation  velocity  due  to  the  moving  frame 
of  reference  given  by 


*  ,  DV  D'Vr 

at  —  {0‘txiatyiatz)  —  jyy 

fliQ 

=  a0  +  (—  X  +  2(0  x  K)  +  0  x  (O  x  r~) 


(15) 


The  contravariant  velocity  components  used  in  Eqs.  (4)  -(7)  are  given  by 


Ur  —  d"  £yVr  d" 

(16) 

Vr  =  V xUr  d-  T]yVr  -f  T)zWr 

(U) 

Wr  —  d"  CyVr  d"  CzV^r 

(18) 

Equation  (2)  is  the  most  general  form  of  thin-layer  NS  equations  in  the  blade-fixed  moving 
frame  of  reference  in  terms  of  relative  velocity.  The  5-term  in  Eq.  (2)  is  the  source  term 

contributed  from  the  time-dependent  rotation  (ft(t))  and  translation  (V0(tj)  of  the  rotor  blade. 

^  -4 

For  the  hovering  motion  case  considered  here,  V0  =  0  and  Q  =  Oi  —  Qj  +  Ok  with  Q,  =  constant, 
the  5-term  in  Eq.  (2),  or  Eq.  (8),  reduces  to 

5  =  p[0,  (x  —  xp)Sl2  +  2  £lwr,  0,  ( z  —  zp)Cl2  —  2fhzr,  0]*  (19) 

For  fixed-wing  motion,  the  source  term  S  is  zero  and  the  relative  velocity  Vr  is  the  same  as  the 
absolute  velocity  V;  in  this  case,  the  TLNS3DR  code  is  identical  to  the  fixed-wing  TLNS3D 
code. 
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2.2  Finite- Volume  Discretization 


The  TLNS3DR  code  is  based  on  the  fixed-wing  TLNS3D  code3  as  mentioned  earlier.  The 
difference  of  the  formulation  of  TLNS3DR  from  that  of  TLNS3D  is  due  to  the  source  term  5 
Hence  the  numerical  scheme  follows  exactly  the  one  of  the  TLNS3D  code  with  the  addition 
ot  the  source  term  5,  along  with  the  necessary  modification  of  the  boundary  conditions  for 
hovering  motion. 


The  semi-discrete  form  of  Eq. 
central  differences, 


(2)  is  obtained  after  replacing  the  spatial  derivatives  by 


fat  (J  Ur)i,j,k  +  (Fri+1/2Jtk  -  Fri-l/2,j,k) 

+  (Gri,j+l/2,k  ~  GriJ- 1/2, k) 

+  {Hri,j,k- 1-1/2  —  Hri'j'k- 1/2  ) 

=  (Gvri,j+l/2,k  ~  Gvri, j-l/2,k)  + 


(20) 


Using  the  notation  Qr  to  represent  all  inviscid  flux  terms,  Qvr  to  represent  the  viscous 
tlux  terms  and  Dr  to  represent  the  artificial  dissipative  fluxes  for  convenience,  Eq  (20)  is 
rewritten  as  ’  v  J 


Gr)i,j,k  +  Qr  —  Qvr  +  Dr  +  Jijk^i,j,k 


(21) 


2.3  Multi-stage  Time-Stepping  Scheme 

rn T  sovlFce  rle1rn\  ?  1S  added  into  the  five-stage  Runge-Kutta  time-stepping  scheme  of  the 

lLJNboU  code.  The  time-stepping  scheme  thus  becomes 


C7(°)  =  jj(n) 


vw  =  vm  _  ai  +  Mi o)  _  llDm  _  j- ls] 


J 
At 


um  =  V(  o)  _  „  +  _  TlX><o)  _  72D(,)  _  j-1|SJ 


J 
At 


P<3>  =  P<0)  -  03  ~[QW  +  MW  -  TlD<»)  -  T2c<»  -  13dW  -  J-'s\ 


J 
At 


uw  =  uw  -  +  MW  -  7,d<°>  -  72d<‘>  -  73dW  -  7iDW  -  J-S] 


(22) 


J 
At 


u(r5)  =  up  -  «5—[qP  +  -  llDp  -  12dP  -  lzDp  -  7,Dp  -  75dP 


-J-'S] 

jj(n+l)  _  jj( 5) 


where  the  superscripts  (n)  and  (n  4-  1)  refer  to  time  levels  and  the  superscripts  (0)  through 
(5)  refer  to  multistage  levels  within  each  time  step. 

In  addition  to  the  five-stage  time-stepping  scheme,  the  local  time-stepping  and  residual 
smoothing  techniques  are  applied  to  accelerate  the  convergence  of  the  scheme  to  steady  state 
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for  both  fixed-wing  and  hovering  motions.  The  details  of  the  residual  smoothing  axe  discussed 
in  Ref.  [3]. 

2.4  Multigrid  Acceleration  Technique 

The  multigrid  scheme  of  TLNS3D  is  employed.  The  idea  of  the  multigrid  is  the  use  of 
successively  coarser  grids  to  compute  corrections  to  a  fine  grid  solution.  The  use  of  multigrid 
technique  can  significantly  accelerate  the  convergence  rate.  The  coefficients  of  ft’s  and  7’s  in 
Eq.  (22)  as  well  as  those  for  the  residual  smoothing  are  carefully  selected  to  provide  optimum 
damping  of  errors.  The  results  presented  in  this  report  are  obtained  using  a  V-cycle  multigrid 
scheme. 

2.5  Boundary  Conditions  and  Grid 

A  C-0  type  grid  is  employed.  The  grid  is  generated  using  either  GRIDGEN3D  or  WTCO3 
code.  On  the  blade  surface,  no-slip  and  no-penetration  conditions  are  used  by  setting  relative 

velocity  (Vr)  to  zero.  The  adiabatic  wall  condition  and  zero-normal  pressure  gradient  condition 
at  the  wall  are  also  applied  at  the  blade  surface.  The  farfield  boundaries  are  treated  by  using 
Riemann  invariants  condition.  At  the  inboard  plane,  near  the  rotation  axis  of  the  blade,  the 
condition  of  ur/wr  =  ut/wt  =  —(z  —  zp)/(x  —  xp )  is  used  to  guarantee  that  the  radial  velocity 
is  zero  at  this  boundary  for  hovering  motion;  for  the  fixed-wing  motion,  the  plane-of-symmetry 
condition  is  used. 
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3.  COMPUTATIONAL  RESULTS 


Application  of  the  TLNS3DR  code  is  made  for  realistic  helicopter  rotor  blade  configu¬ 
rations  to  study  the  rotor  design.  For  nonrotating  flows,  calculated  results  of  the  Berp-tip 
the  swept-tip  and  the  equivalent  taper-tip  are  presented  in  the  Subsection  3.1,  along  with 
comparison  of  some  experimental  data.  For  both  rotating  and  nonrotating  flows,  the  similar 
results  are  presented  in  the  Subsection  3.2  for  the  Berp-type  blade,  the  swept-type  blade,  and 
the  equivalent  taper-type  and  rectangle-type  blades  without  twist. 

3.1  Nonrotating  Blade  Tips 

The  planforms  of  the  Berp-tip,  the  swept-tip,  and  the  taper-tip  are  given  in  Figures  1  2 
respectively.  Figures  4,  5  and  6  are  the  views  of  the  twist  and  airfoil  section  geometry. 
The  C-0  grids  are  generated  for  the  Berp-tip  by  GRIDGEN3D  and  for  the  swept-  and  the 
taper-tips  by  WTCO  grid  generators,  where  a  single  isolated  blade-tip  is  considered  in  each 
case.  Figures  7,  8  and  9  are  the  partial  views  of  the  mesh  around  the  Berp-tip,  the  swept-tip 
and  the  taper-tip,  respectively,  where  129  x  65  x  49  mesh  points  are  used  in  the  streamwise, 
normal  and  spanwise  directions,  respectively.  The  grid  of  this  size  is  used  in  most  cases.  Flows 
at  three  angle-of- attack  conditions  are  considered.  Tip  flow  fields  in  terms  of  the  downwash, 
axial  velocity  and  crossflow  velocity  vector  field  are  presented.  The  surface  pressure  coefficients 
and  the  convergence  history  in  terms  of  the  residual  are  also  given. 

3.1.1  Tip-flow  field  at  a  =  15° 


Figures  10,  11  and  12  are  the  calculated  downwash  velocity  contours  for  the  Berp-tip,  the 
swept-tip  and  the  taper-tip,  respectively,  at  angle  of  attack  of  15°  under  an  incompressible 
tree-stream  Mach  number.  Figure  13  gives  the  same  result  for  the  taper-tip  with  a  finer 
grid  where  161  x  97  x  97  mesh  points  are  used  to  study  the  effect  of  the  mesh  refinement 
to  the  solution.  Figures  14,  15  and  16  are  the  experimentally  measured  downwash  velocity 
contours  for  the  Berp-tip,  the  swept-tip  and  the  taper-tip,  respectively.  The  comparison  of 
the  calculated  results  with  the  experimental  data  shows  that  both  the  TLNS3DR  code  and 
experiment  produce  similar  results  in  the  region  near  the  blade  tips,  and  that  the  tip  vortex 
is  diffused  when  it  convects  to  downstream  in  the  calculated  solution.  The  comparison  of  the 
solution  of  Figure  12  with  that  of  Figure  13  shows  that  the  finer  grid  does  give  a  slightly  better 
solution.  ° 

Figures  17,  18  and  19  are  the  calculated  axial  velocity  contours  for  the  Berp-tip,  the 
swept-tip  and  the  taper-tip,  respectively.  Figure  20  gives  the  same  result  for  the  taper-tip 
with  a  finer  grid  where  161  x  97  x  97  mesh  points  are  used.  Figures  21,  22  and  23  are  the 
experimentally  measured  axial  velocity  contours  for  the  Berp-tip,  the  swept-tip  and  the  taper- 
tip,  respectively.  Both  calculated  solutions  and  experimental  data  give  quantitatively  similar 
results.  It  is  seen  that  the  Berp-tip  produces  the  largest  tip- vortex  in  terms  of  the  diameter  of 
the  vortex  cylinder  from  both  calculated  results  and  experiment.  From  the  calculated  results, 
the  taper-tip  has  a  tightest  tip  vortex  structure;  from  the  experiment,  the  swept-tip  has  the 
tightest  tip  vortex  structure. 

Figures  24,  25  and  26  are  the  calculated  crossflow  velocity  vector  fields  for  the  Berp-tip, 
the  swept-tip  and  the  taper-tip,  respectively.  Figure  27  gives  the  same  result  for  the  taper-tip 
with  a  finer  grid  where  161  x  97  x  97  mesh  points  are  used.  Figures  28,  29  and  30  are  the 
experimentally  measured  crossflow  velocity  vector  fields  for  the  Berp-tip,  the  swept-tip  and  the 
taper-tip,  respectively.  Both  calculated  solution  and  experimental  data  show  quantitatively 
similar  tip  vortex  pattern,  and  tell  us  that  the  strength  of  the  tip  vortex  in  terms  of  the 
crossflow  velocity  of  the  Berp-tip  is  slightly  weaker  than  that  of  the  swept-tip.  The  calculated 


6 


solution  of  the  crossflow  field  of  the  taper-tip  is  likely  diffused. 

3.1.2  Tip-flow  field  at  a  =  5°  —  7.5° 

Figures  31,  32  and  33  are  the  calculated  downwash  velocity  contours  for  the  Berp-tip, 
the  swept-tip  and  the  taper-tip,  respectively,  at  angle  of  attack  of  5°  to  7.5°  under  an  in¬ 
compressible  free-stream  Mach  number.  At  this  angle-of-attack  condition,  the  effective  angle 
of  attack  is  near  zero  due  to  the  negative  twist  angle  near  the  tip.  Figures  34,  35  and  36 
are  the  calculated  axial  velocity  contours  for  the  Berp-tip,  the  swept-tip  and  the  taper-tip, 
respectively.  Figures  37,  38  and  39  are  the  calculated  crossflow  velocity  vector  fields.  Figure 
40  shows  the  experimentally  measured  crossflow  velocity  vector  fields  for  these  three  tips  at 
x  =  1.65  chord  length  station  measured  from  the  leading  edge  of  the  root.  The  calculated 
solution  shows  reasonable  results  in  terms  of  the  downwash,  axial  velocity  and  crossflow  ve¬ 
locity,  but  at  the  location  away  from  the  blade,  such  as  at  x  =  1.65  chord  length  station,  the 
tip  vortex  is  diffused  due  to  the  coarse  grid  used  in  that  region. 

3.1.3  Tip-flow  field  at  a  =  0° 

Figures  41,  42  and  43  are  the  calculated  downwash  velocity  contours  for  the  Berp-tip, 
the  swept-tip  and  the  taper-tip,  respectively,  at  angle  of  attack  of  0°  under  an  incompressible 
free-stream  Mach  number.  At  this  angle-of-attack  condition,  the  effective  angle  of  attack  is  a 
negative  value  due  to  the  negative  twist  angle  near  the  tip.  Figure  44  gives  the  same  result  for 
the  taper-tip  with  a  finer  grid  where  161  X  97  x  97  mesh  points  are  used.  Figures  45,  46  and 
47  are  the  calculated  axial  velocity  contours.  Figure  48  gives  the  same  result  for  the  taper-tip 
with  a  finer  grid  where  161  x  97  X  97  mesh  points  are  used.  Figures  49,  50  and  51  are  the 
calculated  crossflow  velocity  vector  fields.  Figure  52  gives  the  same  result  for  the  taper-tip 
with  a  finer  grid  where  161  x  97  x  97  mesh  points  are  used.  Figures  53  is  the  experimentally 
measured  crossflow  velocity  vector  fields  for  the  Berp-tip,  the  swept-tip  and  the  taper-tip  at 
x-station  of  1.65  chord  length  from  the  leading  edge  of  the  root.  Calculated  solution  and 
experimental  data  agree  quantitatively.  It  can  be  seen  from  both  the  experimental  data  and 
the  calculated  solution  that  the  Berp-tip  does  not  produce  a  significant  tip  vortex,  and  that 
the  tip  vortex  from  the  taper-tip  is  flatter  than  that  from  the  swept-tip.  Since  the  effective 
angle  of  attack  is  a  negative  value,  the  tip  vortex  rotates  in  the  clockwise  direction  and  formed 
under  the  blade  tip. 

3.1.4  Surface  pressure  and  convergence  history 

Figures  54,  55  and  56  are  the  calculated  surface  pressure  coefficient  contours  for  the  Berp- 
tip,  the  swept-tip  and  the  taper-tip,  respectively,  at  angle  of  attack  of  15°.  Figures  57,  58 
and  59  are  the  calculated  surface  pressure  coefficient  contours  at  angle  of  attack  of  5°  to  7.5°. 
Figures  60,  61  and  62  are  the  calculated  surface  pressure  coefficient  contours  at  angle  of  attack 
of  0°.  These  results  are  self-explanatory. 

Figures  63,  64  and  65  are  the  convergence  history  in  terms  of  residual  for  the  Berp-tip, 
the  swept-tip  and  the  taper-tip,  respectively.  Figure  66  is  the  convergence  history  for  the 
taper-tip  with  a  finer  grid.  As  mentioned  earlier,  a  V-cycle  multigrid  scheme  is  employed.  A 
total  of  375  equivalent  fine  grid  time  steps  (based  on  the  equivalent  CPU  time  required)  are 
used.  The  CPU  time  on  the  NASA-Langley’s  Cray-YMP  with  a  single  processor  for  cases 
with  129  x  65  x  49  fine  grid  points  is  4600  seconds,  which  is  0.011  second  per  fine  grid  point, 
or  29.8  x  10-6  second  per  fine  grid  per  time  step.  Comparing  with  the  single-grid  method,  the 
current  code  is  very  efficient.  Figures  65  and  66  show  that,  on  the  other  hand,  the  convergence 
rate  is  dependent  on  the  grid  size,  and  that  the  rate  decreases  when  the  grid  size  decreases. 
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3.2  Comparison  of  Rotating  and  Nonrotating  Rotor  Flows 

The  planforms  and  partial  views  of  the  C-0  grid  of  the  Berp-blade,  the  swept-blade,  and 
the  equivalent  taper-blade  and  rectangle-blade  without  twist  are  given  in  Figures  67  through 
74.  A  mesh  of  129  x  65  x  65  grid  points  is  used  in  the  streamwise,  normal  and  spanwise 
directions,  respectively.  Flows  at  angle  of  attack  of  2  degrees  under  a  compressible  free- 
stream  /  tip  Mach  number  of  0.628  are  considered  for  both  nonrotating  and  rotating  flows. 
Tip  flow  fields  in  terms  of  the  axial  velocity  and  crossflow  velocity  vector  field  are  presented. 
The  blade-surface  pressure  coefficients,  local  Mach  number,  local  spanwise  velocity  and  the 
convergence  history  in  terms  of  the  residual  are  also  presented. 

For  nonrotating  motion  of  the  Berp-blade,  calculated  tip  flow  field  in  terms  of  axial 
velocity  and  crossflow  velocity  vector  are  presented  in  Figures  75  and  76,  respectively.  Figures 
77,  78  and  79  are  the  calculated  surface  pressure  coefficients,  surface  local  Mach  number  and 
local  spanwise  velocity,  respectively,  where  the  surface  local  Mach  number  and  local  spanwise 
velocity  are  calculated  at  j  =  33  which  is  about  just  outside  of  the  boundary  layer.  Figure  80 
is  the  convergence  history. 

For  rotating,  hovering  motion  of  the  Berp-blade,  calculated  tip  flow  field  in  terms  of 
relative  and  absolute  axial  velocities  and  relative  and  absolute  crossflow  velocity  vectors  are 
presented  in  Figures  81,  82,  83  and  84,  respectively.  Figures  85,  86,  87  and  88  are  calculated 
surface  pressure  coefficients,  surface  local  relative  Mach  number  and  local  relative  and  absolute 
spanwise  velocities,  respectively.  Figure  89  is  the  convergence  history. 

The  same  sets  of  results  as  those  for  the  Berp-blade  are  presented  for  the  swept-blade,  the 
taper-blade  and  the  rectangle-blade.  Figures  90  through  104  are  results  of  the  swept-blade; 
Figures  105  through  119  are  results  of  the  taper-blade;  and  Figures  120  through  134  are  results 
of  the  rectangle-blade. 

All  results  presented  are  self-explanatory.  The  CPU  time  used  is  6100  seconds  on  the 
Cray-YMP  with  a  single  processor.  The  results  on  the  rotating,  hovering  motion  of  the  blade 
are  quantitatively  correct  where  tip  vortex  in  terms  of  both  relative  and  absolute  velocities  are 
captured,  and  where  the  blade  surface  flow  field  shows  that  the  rotating  effects  are  correctly 
calculated. 
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4.  CONCLUDING  REMARKS 


The  TLNS3DR  code  is  applied  to  various  realistic  helicopter  rotor  blade  configurations 
to  study  the  advanced  rotor  tip  flow  and  to  further  investigate  the  capability  of  the  code.  For 
nonrotating  tip  geometry,  three  tips  are  computationally  investigated  and  briefly  compared 
with  the  experimental  data.  For  both  rotating  hovering  motion  and  nonrotating  motion  of 
blade  geometry,  four  blades  are  computationally  investigated.  The  results  of  rotating  flows  are 
presented  along  with  those  of  equivalent  nonrotating  results  for  studying  the  effects  of  rotation 
term  in  the  TLNS3DR  code  to  validate  the  code.  The  detailed  results  with  limited  discussion 
presented  in  this  report  are  expected  to  be  useful  in  undestanding  the  rotor  design,  particularly 
in  studying  the  tip  vortex.  Througth  this  investigation,  several  concluding  remarks  can  be 
drawn: 

1.  The  TLNS3DR  code  is  an  efficient  and  accurate  CFD  tool  for  nonrotating  flow  calculations. 

2.  The  TLNS3DR  code  efficiently  produces  quantitatively  correct  solution  for  rotating  hovering 
motion;  however,  the  accuracy  of  the  code  for  rotating  flows  needs  to  be  tested. 

3.  The  TLNS3DR  code  has  capability  to  capture  the  tip  vortex  and  the  solution  of  the  tip 
vortex  generally  agrees  with  that  of  experiment;  it  is  very  promising  for  using  the  code  as  an 
analysis  tool  in  rotor  tip  design. 

4.  The  structure  of  the  tip  vortex  does  depend  on  the  planform  of  the  rotor  blade,  and  it  is 
predictable  by  the  TLNS3DR  code. 

5.  A  computational  grid  finer  than  the  current  ones  is  needed  to  accurately  capture  the  tip 
vortex. 
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Plane  Form  of  BART  Berp-type  Wing/Rotor  Tip,  1 29x65x49  Mesh 


11 


Plane  Form  of  BART  Swept-type  Wing/Rotor  Tip,  1 29x65x49  Mesh 

Figure  2.  Plane  form  of  swept-tip. 
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Plane  Form  of  BART  Taper-type  Wing/Rotor  Tip,  1 29x65x49  Mesh 

Figure  3.  Plane  form  of  taper-tip. 
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Z-coordinate  is  expanded  by  14  times. 

View  of  The  Twist  for  BART  Berp-type  Wing/Rotor  Tip,  1 29x65x49  Mesh 
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Z-coordinate  is  expanded  by  14  times. 

View  of  The  Twist  for  BART  Swept-type  Wing/Rotor  Tip,  1 29x65x49  Mesh 

Figure  5.  View  of  the  twist  for  swept-tip. 
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Z-coordinate  is  expanded  by  1 4  times. 

View  of  The  Twist  for  BART  Taper-type  Wing/Rotor  Tip,  1 29x65x49  Mesh 

Figure  6.  View  of  the  twist  for  taper-tip. 
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Partial  View  of  C-0  Mesh  for  BART  Berp-type  Wing/Rotor  Tip,  129x65x49  Mesh 
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Partial  View  of  C-0  Mesh  for  BART  Swept-type  Wing/Rotor  Tip,  129x65x49  Mesh 
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Partial  View  of  C-0  Mesh  for  BART  Taper-type  Wing/Rotor  Tip,  1 29x65x49  Mesh 
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BART  Swept-type  Wing  Tip,  a  =15°,  Re=2x106, 129x65x49  Mesh 
Tip  Flow  Field  Investigation  -  Calculated  Downwash  (W)  Contours 

Figure  11.  Calculated  tip  flow  field  in  terms  of  downwash  for  swept-tip  at  a=15°. 
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BART  Taper-type  Wing  Tip,  a  =  1 5°,  Re=2x1 06, 161  x97x97  Mesh 
Tip  Flow  Field  Investigation  -  Calculated  Downwash  (W)  Contours 

Figure  13.  Calculated  tip  flow  field  in  terms  of  downwash  for  taper-tip  at  a=15°  with  a  fine  grid. 
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Figure  15.  Measured  tip  flow  field  in  terms  of  downwash  for  swept-tip  at  a=15°. 
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Swept-type  Wing  Tip,  a  =1 5°,  Re=2x1 06, 1 29x65x49  Mesh 
Field  Investigation  -  Calculated  Axial  (U-)  Velocity  Contours 

Figure  18.  Calculated  tip  flow  field  in  terms  of  axial  velocity  for  swept-tip  at  a=15°. 


Figure  20.  Calculated  tip  flow  field  in  terms  of  axial  velocity  f< 
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srp-type  Wing  Tip,  a  =15°,  Re=2x106, 129x65x49  Mesh 
'  Field  Investigation  -  Calculated  V-W  Velocity  Vectors 

Calculated  tip  flow  field  in  terms  of  crossflow  velocity  vector  for  Berp-tip  at  a=15°. 
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BART  Swept-type  Wing  Tip,  a  =15°,  Re=2x106, 129x65x49  Mesh 
Tip  Flow  Field  Investigation  -  Calculated  V-W  Velocity  Vectors 

Figure  25.  Calculated  tip  flow  field  in  terms  of  crossflow  velocity  vector  for  swept-tip  at  a=15°. 
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BART  Taper-type  Wing  Tip,  a  =15°,  Re=2x106, 129x65x49  Mesh 
Tip  Flow  Field  Investigation  -  Calculated  V-W  Velocity  Vectors 

Figure  26.  Calculated  tip  flow  field  in  terms  of  crossflow  velocity  vector  for  taper-tip  at  a=15°. 
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BART  Taper-type  Wing  Tip,  a  =15°,  Re=2x106, 161x97x97  Mesh 
Tip  Flow  Field  Investigation  -  Calculated  V-W  Velocity  Vectors 

Figure  27.  Calculated  tip  flow  field  in  terms  of  crossflow  velocity  vector  for  taper-tip  at  a=15°  with  a  fine  grid. 
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Figure  29.  Measured  tip  flow  field  in  terms  of  crossflow  velocity  vector  for  swept-tip  at  a=15°. 
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BART  Berp-type  Wing  Tip,  a  =5°,  Re=2x1 06, 1 29x65x49  Mesh 
Tip  Flow  Field  Investigation  -  Calculated  Downwash  (W)  Contours 

Figure  31.  Calculated  tip  flow  field  in  terms  of  downwash  for  Berp-tip  at  a=5°. 
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•type  Wing  Tip,  a  =5°,  Re=2x1 06, 1 29x65x49  Mesh 
Investigation  -  Calculated  Downwash  (W)  Contours 

Calculated  tip  flow  field  in  terms  of  downwash  for  swept-tip  at  a=5°. 
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BART  Taper-type  Wing  Tip,  a  =7.5°,  Re=2x106, 129x65x49  Mesh 
Tip  Flow  Field  Investigation  -  Calculated  Downwash  (W)  Contours 

Figure  33.  Calculated  tip  flow  field  in  terms  of  downwash  for  taper-tip  at  a=7.5°. 
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BART  Swept-type  Wing  Tip,  a  =5°,  Re=2x1 06, 1 29x65x49  Mesh 
Tip  Flow  Field  Investigation  -  Calculated  Axial  (U-)  Velocity  Contours 

Figure  35.  Calculated  tip  flow  field  in  terms  of  axial  velocity  for  swept-tip  at  a=5°. 


ir-type  Wing  Tip,  a  =7.5  ,  Re=2x10 , 129x65x49  Mesh 
d  Investigation  -  Calculated  Axial  (U-)  Velocity  Contours 

Calculated  tip  flow  field  in  terms  of  axial  velocity  for  taper-tip  at  a=7.5°. 


47 


BART  Swept-type  Wing  Tip,  a  =5°,  Re=2x1 06, 1 29x65x49  Mesh 
Tip  Flow  Field  Investigation  -  Calculated  V-W  Velocity  Vectors 

Figure  38.  Calculated  tip  flow  field  in  terms  of  crossflow  velocity  vector  for  swept-tip  at  a=5°. 
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BART  Taper-type  Wing  Tip,  a  =7.5°,  Re=2x106, 129x65x49  Mesh 
Tip  Flow  Field  Investigation  -  Calculated  V-W  Velocity  Vectors 

Figure  39.  Calculated  tip  flow  field  in  terms  of  crossflow  velocity  vector  for  taper-tip  at  a=7.5°. 


50 


BART  Berp-type  Wing  Tip,  a  =0°,  Re=2x1 06, 1 29x65x49  Mesh 
Tip  Flow  Field  Investigation  -  Calculated  Downwash  (W)  Contours 

Figure  41.  Calculated  tip  flow  field  in  terms  of  downwash  for  Berp-tip  at  a=0°. 
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BART  Taper-type  Wing  Tip,  a  =0°,  Re=2x1 06, 1 29x65x49  Mesh 
Tip  Flow  Field  Investigation  -  Calculated  Downwash  (W)  Contours 

Figure  43.  Calculated  tip  flow  field  in  terms  of  downwash  for  taper-tip  at  a=0°. 
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BART  Taper-type  Wing  Tip,  a  =0°,  Re=2x106,  161x97x97  Mesh 
Tip  Flow  Field  Investigation  -  Calculated  Downwash  (W)  Contours 

Figure  44.  Calculated  tip  flow  field  in  terms  of  downwash  for  taper-tip  at  a=0°  with  a  fine  grid. 
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BART  Berp-type  Wing  Tip,  a  =0°,  Re=2x106, 129x65x49  Mesh 
Tip  Flow  Field  Investigation  -  Calculated  Axial  (U-)  Velocity  Contours 

Figure  45.  Calculated  tip  flow  field  in  terms  of  axial  velocity  for  Berp-tip  at  a=0°. 
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BART  Swept-type  Wing  Tip,  a  =0°,  Re=2x1 06, 1 29x65x49  Mesh 
Tip  Flow  Field  Investigation  -  Calculated  Axial  (U-)  Velocity  Contours 

Figure  46.  Calculated  tip  flow  field  in  terms  of  axial  velocity  for  swept-tip  at  cc=0°. 
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BART  Taper-type  Wing  Tip,  a  =0°,  Re=2x106, 1 29x65x49  Mesh 
Tip  Flow  Field  Investigation  -  Calculated  Axial  (U-)  Velocity  Contours 

Figure  47.  Calculated  tip  flow  field  in  terms  of  axial  velocity  for  taper-tip  at  a=0°. 
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BART  Taper-type  Wing  Tip,  a  =0°,  Re=2x106, 161x97x97  Mesh 
Tip  Flow  Field  Investigation  -  Calculated  Axial  (U-)  Velocity  Contours 

Figure  48.  Calculated  tip  flow  field  in  terms  of  axial  velocity  for  taper-tip  at  a=0°  with  a  fine  grid. 
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BART  Berp-type  Wing  Tip,  a  =0°,  Re=2x106, 129x65x49  Mesh 
Tip  Flow  Field  Investigation  -  Calculated  V-W  Velocity  Vectors 

Figure  49.  Calculated  tip  flow  field  in  terms  of  crossflow  velocity  vector  for  Berp-tip  at  a=0°. 
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BART  Swept-type  Wing  Tip,  a  =0°,  Re=2x1 06, 1 29x65x49  Mesh 
Tip  Flow  Field  Investigation  -  Calculated  V-W  Velocity  Vectors 

Figure  50.  Calculated  tip  flow  field  in  terms  of  crossflow  velocity  vector  for  swept-tip  at  a=0°. 
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BART  Taper-type  Wing  Tip,  a  =0°,  Re=2x106, 161x97x97  Mesh 
Tip  Flow  Field  Investigation  -  Calculated  V-W  Velocity  Vectors 

Figure  52.  Calculated  tip  flow  field  in  terms  of  crossflow  velocity  vector  for  taper-tip  at  a=0°  with  a  fine  grid. 


BART  BERP-type  WING  TIP.  a=0,  X/ C=1.65  BART  SWEPT-type  TIP.  a=0,  X/C=1.65 


o 


62 


Figure  53.  Measured  tip  flow  field  in  terms  of  crossflow  velocity  vector  for  three  tips  at  a=0°  and  x/c=1.65. 
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Berp-type  Wing  Tip,  ct=15°,  Re=2x106, 129x65x49  Mesh 
Calculated  Surface  Pressure  Coefficient  Contours 

Figure  54.  Calculated  surface  pressure  coefficients  for  Berp-tip  at  a=15°. 
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BART  Swept-type  Wing  Tip,  ct=  1 5°,  Re=2x1 0®,  1 29x65x49  Mesh 
Calculated  Surface  Pressure  Coefficient  Contours 

Figure  55.  Calculated  surface  pressure  coefficients  for  swept-tip  at  a=15°. 
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BART  Taper-type  Wing  Tip,  <x=15°,  Re=2x106, 129x65x49  Mesh 
Calculated  Surface  Pressure  Coefficient  Contours 

Figure  56.  Calculated  surface  pressure  coefficients  for  taper-tip  at  a=15°. 
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BART  Berp-type  Wing  Tip,  a=5°,  Re=2x106, 129x65x49  Mesh 
Calculated  Surface  Pressure  Coefficient  Contours 
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BART  Swept-type  Wing  Tip,  a=5°,  Re=2x1 06, 1 29x65x49  Mesh 
Calculated  Surface  Pressure  Coefficient  Contours 

Figure  58.  Calculated  surface  pressure  coefficients  for  swept-tip  at  a=5°. 
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BART  Taper-type  Wing  Tip,  a=7.5°,  Re=2x1 06, 1 29x65x49  Mesh 
Calculated  Surface  Pressure  Coefficient  Contours 

Figure  59.  Calculated  surface  pressure  coefficients  for  taper-tip  at  a=7.5°. 
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Berp-type  Wing  Tip,  <x=0°,  Re=2x1 06, 1 29x65x49  Mesh 
Calculated  Surface  Pressure  Coefficient  Contours 

Figure  60.  Calculated  surface  pressure  coefficients  for  Berp-tip  at  a=0°. 
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BART  Swept-type  Wing  Tip,  cx=0o,  Re=2x1 06, 1 29x65x49  Mesh 
Calculated  Surface  Pressure  Coefficient  Contours 

Figure  61.  Calculated  surface  pressure  coefficients  for  swept-tip  at  a=0°. 
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BART  Taper-type  Wing  Tip,  a=0°,  Re=2x1 06, 1 29x65x49  Mesh 
Calculated  Surface  Pressure  Coefficient  Contours 

Figure  62.  Calculated  surface  pressure  coefficients  for  taper-tip  at  a=0°. 
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Figure  64.  Convergence  history  in  terms  of  residuals  for  swept-tip  at  a=15°,  a=5°  and  a=0°. 
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Plane  Form  of  Berp-type  Blade  with  0=2°,  129x65x65  Mesh 
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Plane  Form  of  Swept-type  Blade  with  0=2°,  1 29x65x65  Mesh 
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Partial  View  of  C-0  Mesh  for  Swept-type  Blade  with  9=2°,  129x65x65  Mesh 
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Plane  Form  of  Taper-type  Blade  with  0=2°,  129x65x65  Mesh 


81 


Partial  View  of  C-0  Mesh  for  Taper-type  Blade  with  0=2°,  129x65x65  Mesh 


82 


Plane  Form  of  Rectangle-type  Blade  with  0=2°,  1 29x65x65  Mesh 

Figure  73.  Plane  form  of  rectangle-blade. 
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Partial  View  of  C-0  Mesh  for  Rectangle-type  Blade  with  6=2°,  1 29x65x65  Mesh 


Nonrotating  Berp-type  Blade 
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e=a  =2°,  Mx=. 628,  Re=3.7x106, 129x65x65  Mesh 
Tip  Flow  Field  Investigation  -  Calculated  Axial  (U-)  Velocity  Contours 

Figure  75.  Calculated  tip  flow  field  in  terms  of  axial  velocity  for  nonrotating  Berp-blade. 
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6=a=2°,  M„=.628,  Re=3.7x106, 129x65x65  Mesh 
Calculated  Surface  Pressure  Coefficient  Contours 
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M^=.628,  Re=3.7x1 06, 129x65x65  Mesh 
Calculated  Local  Mach  Contours 


Figure  79.  Calculated  local  V-velocity  contours  for  nonrotating  Berp-blade. 


(lenpisay)  0l6cr| 


CO 

C/3 

D 

0 

T3 

W 

0 

in 

CD 

DC 

X 

in 

M— 

o 

CD 

C/3 

X 

E 

03 

1— 

CM 

CD 

T- 

c 

CD 

■ 

o 

>■, 

y— 

l— 

X 

o 

N- 

w 

00 

II 

X 

1 1 

0 

DC 

O 

c 

co~ 

0 

CM 

03 

CD 

k_ 

0 

ll" 

> 

8 

C 

o 

O 

89 


Figure  80.  Convergence  history  in  terms  of  residuals  for  nonrotating  Berp-blade. 


Berp-type  Blade  in  Hover 


as 


0=a  =2  ,  Mtip=.628,  Q=0.0551 ,  Re=3.7x1 06,  65x33x33  Mesh 
Tip  Flow  Field  Investigation  -  Calculated  Relative  Axial  Velocity  Contours 

Figure  81.  Calculated  tip  flow  field  in  terms  of  relative  axial  velocity  for  Berp-blade  in  hover. 


Berp-type  Blade  in  Hover 


Mtlp=.628,  Q=0.0551 ,  Re=3.7x106,  65x33x33  Mesh 
Investigation  -  Calculated  Absolute  Axial  Velocity  Contours 

Calculated  tip  flow  field  in  terms  of  absolute  axial  velocity  for  Berp-blade  in  hover. 
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Berp-type  Blade  in  Hover 


93 


(Interpolated  to  1 29x65x65  Mesh)  ^w^i.es  x/c=i.25 

e=«  .=2°.M,iP=-628,  Q=0.0551 ,  Re=3.7x106,  65x33x33  Mesh 
Tip  Flow  Field  Investigation  -  Calculated  Absolute  V-W  Velocity  Vectors 

Figure  84.  Calculated  tip  flow  field  in  terms  of  absolute  crossflow  velocity  vector  for  Berp-blade  in  hover. 
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=2  ,  Mtlp=.628,  Q=0.0551 ,  Re=3.7x106,  65x33x33  Mesh 
Calculated  Surface  Pressure  Coefficient  Contours 


/Itip=.628,  £2=0.0551,  Re=3.7x106,  65x33x33  Mesh 
Calculated  Local  Relative  Mach  Contours 
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-  MtiP=-628,  0=0.0551 ,  Re=3.7x106,  65x33x33  Mesh 
Calculated  Local  Relative  V-Velocity  Contours 
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Convergence  History  in  terms  of  Residual 


Nonrotating  Swept-type  Blade 
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e=a=2°,  M„=.628,  Re=3.7x106, 129x65x65  Mesh 
Calculated  Surface  Pressure  Coefficient  Contours 

Figure  92,  Calculated  surface  pressure  coefficients  for  nonrotating  swept-blade. 
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M„=.628,  Re=3.7x1 06, 129x65x65  Mesh 
Calculated  Local  Mach  Contours 
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M„=.628,  Re=3.7x1 06, 129x65x65  Mesh 
Calculated  V-Velocity  Contours 
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Figure  95.  Convergence  history  in  terms  of  residuals  for  nonrotating  swept-blade. 


Swept-type  Blade  in  Hover 


105 


\1X/C=1.6S,a'v=' 

0=a  =2°,  Mtlp=.628,  Q=0.0479,  Re=3.7x106, 129x65x65  Mesh 
Tip  Flow  Field  Investigation  -  Calculated  Relative  Axial  Velocity  Contours 

Figure  96.  Calculated  tip  flow  field  in  terms  of  relative  axial  velocity  for  swept-blade  in  hover. 
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e=a  =2  ,  Mtip=.628,  Q=0.0479,  Re=3.7x1 06, 1 29x65x65  Mesh 
Tip  Flow  Field  Investigation  -  Calculated  Relative  V-W  Velocity  Vectors 

Figure  98.  Calculated  tip  flow  field  in  terms  of  relative  crossflow  velocity  vector  for  swept-blade  in  hover. 


Swept-type  Blade  in  Hover 
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xjx/c=1 ,65^</C=1 .25 - ~ 

=2°>  MtiP=-628,  Q=0.0479,  Re=3.7x1 06,  1 29x65x65  Mesh 
Tip  Flow  Field  Investigation  -  Calculated  Absolute  V-W  Velocity  Vectors 

Figure  99.  Calculated  tip  flow  field  in  terms  of  absolute  crossflow  velocity  vector  for  swept-blade  in  hover. 
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0=a=2  ,  Mt|p=.628,  Q=0.0479,  Re=3.7x106,  129x65x65  Mesh 
Calculated  Local  Relative  Mach  Contours 


atj=33 
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0=a=2°,  Mtip=.628,  Q=0.0479,  Re=3.7x106, 129x65x65  Mesh 
Calculated  Local  Relative  V-Velocity  Contours 


at  1=33 
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8  a  2  ,  Mljp-.628,  -2-0.0479,  Re=3.7x106, 129x65x65  Mesh 
Calculated  Local  Absolute  V-Velocity  Contours 
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Figure  104.  Convergence  history  in  terms  of  residuals  for  swept-blade  in  hover. 


Nonrotating  Taper-type  Blade 
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X-  rr.  f  ■  M»=-628>  Re=3.7x106, 129X65X65  Mesh 

Tip  Flow  Field  Investigation  -  Calculated  Axial  (U-)  Velocity  Contours 

Figure  105.  Calculated  tip  flow  field  in  terms  of  axial  velocity  for  nonrotating  taper-blade. 


Nonrotating  Taper-type  Blade 
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e=a  =2  ,  Moo=.628,  Re=3.7x106, 129x65x65  Mesh 
Tip  Flow  Field  Investigation  -  Calculated  V-W  Velocity  Vectors 


Level  NCP 
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0=a=2  ,  M„=.628,  Re=3.7x106, 129x65x65  Mesh 
Calculated  Surface  Pressure  Coefficient  Contours 

Figure  107.  Calculated  surface  pressure  coefficients  for  nonrotating  taper-blade. 
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Figure  108.  Calculated  local  Mach  contours  for  nonrotating  taper-blade. 
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0=a=2°,  M„=.628,  Re=3.7x106, 129x65x65  Mesh 
Calculated  V-Velocity  Contours 
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Figure  110.  Convergence  history  in  terms  of  residuals  for  nonrotating  taper- 


Taper-type  Blade  in  Hover 
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e=a  =2  ,  Mtlp=.628,  Q=0.0712,  Re=3.7x106, 129x65x65  Mesh 
Tip  Flow  Field  Investigation  -  Calculated  Relative  Axial  Velocity  Contours 

Figure  111.  Calculated  tip  flow  field  in  terms  of  relative  axial  velocity  for  taper-blade  in  hover. 
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Calculated  tip  flow  field  in  terms  of  absolute  axial  velocity  for  taper-blade  in  hover. 


Taper-type  Blade  in  Hover 


122 


e=a  =2°,  Mtip=.628,  Q=0.0712,  Re=3.7x106, 129x65x65  Mesh 
Tip  Flow  Field  Investigation  -  Calculated  Relative  V-W  Velocity  Vectors 

Figure  113.  Calculated  tip  flow  field  in  terms  of  relative  crossflow  velocity  vector  for  taper-blade  in  hover. 
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Level  NCP 
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2°,  Mtip=.628,  Q=0.071 2,  Re=3.7x106, 129x65x65  Mesh 
Calculated  Surface  Pressure  Coefficient  Contours 
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Level  V 
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M,ip=.628,  Q=0.0712,  Re=3.7x106, 129x65x65  Mesh 
Calculated  Local  Absolute  V-Velocity  Contours 
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Figure  119.  Convergence  history  in  terms  of  residuals  for  taper-blade  in  hover. 


Nonrotating  Rectangle-type  Blade 
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Nonrotating  Rectangle-type  Blade 
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0=a=2°,  Mx=.628,  Re=3.7x106, 129x65x65  Mesh 
Calculated  Surface  Pressure  Coefficient  Contours 

Figure  122.  Calculated  surface  pressure  coefficients  for  nonrotating  rectangle-blade. 


Level  Mach 


<D 

3 

-SP 

E 


132 


M^.628,  Re=3.7x106, 129x65x65  Mesh 
Calculated  Local  Mach  Contours 
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Figure  124.  Calculated  local  V-velocity  contours  for  nonrotating  rectangle-blade. 
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Convergence  History  in  terms  of  Residual 


Rectangle-type  Blade  in  Hover 
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Rectangle-type  Blade  in  Hover 
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^X/C=1.65  X/C=125 

0=a  =2°,  Mtip=.628,  Q=0.0589,  Re=3.7x106, 129x65x65  Mesh 
Tip  Flow  Field  Investigation  -  Calculated  Relative  V-W  Velocity  Vectors 

Figure  128.  Calculated  tip  flow  field  in  terms  of  relative  crossflow  velocity  vector  for  rectangle-blade  in  hover. 


Rectangle-type  Blade  In  Hover 
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^*X/C=1 .65  X/C=1  25 

e=a  =2°,  Mtip=.628,  Q=0.0589,  Re=3.7x106, 129x65x65  Mesh 
Tip  Flow  Field  Investigation  -  Calculated  Absolute  V-W  Velocity  Vectors 

Figure  129.  Calculated  tip  flow  field  in  terms  of  absolute  crossflow  velocity  vector  for  rectangle-blade  in  hover. 
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Level  Mach 
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at  j=33 
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Calculated  Local  Relative  V-Velocity  Contours 
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M^.628,  Q=0.0589,  Re=3.7x106, 129x65x65  Mesh 
Calculated  Local  Absolute  V-Velocity  Contours 
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Figure  134.  Convergence  history  in  terms  of  residuals  for  rectangle-blade  in  hover. 
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SYMBOLS 


at 

Cp 

c 

Dr 

6r 

Fr 

Gr 

Gvr 

Hr 

h  r 

J 

FI  tip 

Moc 

p 

Qr 
Q  vr 

Re 

r  =  (x,y,z) 
rp  =  ( xpi  Vp,  zp) 

s 

t 

Ur 

ur ,  vr:  wr 
V  =  (u,  v ,  w ) 

Ye  =  ( ue,ve,we ) 
K  =  (u0,t;0,u;0) 


acceleration  of  transformation  velocity 
surface  pressure  coefficient 
blade  root  chord  length 
artificial  dissipative  fluxes 
total  energy 

(-component  inviscid  flux 
//-component  inviscid  flux 
//-component  viscous  flux 
(-component  inviscid  flux 
total  enthalpy 
transformation  Jacobian 
tip  Mach  number 
free-stream  Mach  number 
pressure 

all  inviscid  flux  terms 
all  viscous  flux  terms 
Reynolds  number 
position  vector 

pivot  point  vector  or  blade  axis 
source  term  due  to  moving  frame 
time 

field  vector 

contravariant  velocity  components 
absolute  velocity 
rotation  velocity 
translation  velocity 
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VT  —  ( ur,vr,wT )  relative  velocity 

Vt  =  (ut,vt,wt)  transformation  velocity 


a 

7 

e 

e 

B 

Hi  Vi  C) 

p 

Q, 


angle  of  attack 

gas  specific  heat  ratio 

turbulent  eddy-viscosity  correction 

turbulent  eddy-conductivity  correction 

same  as  a 

laminar  viscous  coefficient 
body-fitted  coordinates 
density 

angular  velocity 
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